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by 

Wayne J. Book 


Abstract 


In the second quarterly report of Contract NAS8-28055 a mathematical 
procedure using 4x4 transformation matrices for analyzing the vibrations 
of flexible manipulators was reported and applied to a specific exan^le. 
This report summarizes the previous work and extends the method to include 
flexible Joints as well as links, and to account for the effects of various 
power transmission schemes. A systematic study of the allocation of 
structural material and the placement of components such as motors and 
gearboxes has been undertaken using the tools developed. As one step in 
this direction the variables which relate the vibration parameters of the 
arm to the task and environment of the arm have been Isolated and non- 
dlmenslonallzed. In this manner one is able to reduce the number of 
variables and yet hepefully retain an Intuitive feel for the problem. 

This effort is being continued as a general problem, making reasonable 
assumptions as to the configuration and parameter ranges which are of 
interest to further reduce the large number of variables and arrive at a 
meaningful design tradeoff study. It is desirable at some future point to 
consider a more specific case, whether this case is established by NASA or 
assumed by the investigators. 

The 4x4 transformation matrices have also been used to develop 
analytical expressions for the terms of the complete 6x6 compliance matrix 
for the case of two flexible links joined by a rotating joint, flexible 
about its axis of rotation. This seems to be the most frequently re- 
curring configuration. The availability of these analytical expressions 
in terms of the link and joint parameters will circumvent the numerical 
evaluation of these terms in further studies of this case. 
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The Compliance of Jointed Beams -Practical Matrix Approach 

Manipulator arras are subject to deflection under loads and to vibrations 
about an equilibrium position when the loading on the arm is suddenly 
changed. The deflections deteriorate end point accuracy as computed 
from joint positions and the vibrations can seriously deteriorate the 
response of the arm, and the ability of an operator to perform desired 
maneuvers. The following is a method for analyzing the deflection of 
an arm under given loading conditions. The arm compliance matrix is 
arrived at giving three displacements and three rotations as a linear 
function of the applied forces and moments. The method can be used 
to evaluate the bending of the arm segments and flexible joints as 
well. If the compliance matrix is nonsingular it can be inverted to 
yield a spring constant matrix and hence forces end moments as a 
function of displacements and rotations. The motion of a lumped mass 
spring system can be described by a linear differential equation using 
these spring constants. The validity of this approximation for an arm 
vibrating about an equilibrium position depends largely on how well the 
mass Involved can be lumped into a reasonable number of masses. It is 
less seriously limited by a small amplitude assumption, the assumption 
of negligible damping (only second order effects on the natural frequency), 
and the assumptions that the joint angles are not changing, l^fhen the 
mass of the payload is large compared to the mass of the arm the approxi- 
mation is very good. 


The Mechanics of Arm Deflection 


Consider an arm in static equilibrium with the forces and moments 
on its two ends as is shown in Figure 1. Initially we will assume 

1) the mass of the arm can be lumped for purposes of vibration studies 

2) the arm joints are rigid 

3) the arms segments are simple beams 
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l^hen loads are applied to the ends of the arm the individual arm 
segments deform according to the forces and moments placed on them 
by the neighboring segments. When these forces are expressed in 
terms of a coordinate system which has one axis coincident with 
the neutral axis of the beam as shown in Fig. 2, the deflections 
over the length of the segment are simply obtained . Each of the 
deflections and angles along the three mutually perpendicular 
directions is a linear function of at most two of the loads. 

Notice that one end of the beam is assuraed at the zero position: 


Eq, 1 


Ax 


Ay 

AZ 


e 


X 



0 


z 


“zF “zM 

= 

^z “e^M ^ 


For a beam whose cross section is symmetric about the Y and Z axes 
= oi__ and will be denoted This will be the only case 

specifically considered. The development which follows could retain 
the extra subscript at some loss in readability. The simplification in 
notation is as follows: 


“XF “ “yF “ 

“zF 

“XM “ “yM “ 

®ZM 

“of “ “e^F “ 

“0/ 

“0M “ “e^M “ 
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Beam theory additionally requires that = ^‘qF* ^determining the 

end displacement is a matter of summing the displacements of the in- 
dividual segments and in accounting for the displacement due to end 
point rotations at a distance from the end of the segment where 
the rotation is calculated. For numerous arbitrary joint angles 
this becomes a complex bookkeeping task. The matrix procedure which 
is developed here automatically performs this task. 


Transformation of Coordinates Using 4x4 Matrices 


We are interested in a transformation between two coordinate 
systems whose origins are displaced from one another and whose axes are 
not parallel, as in Fig. 3. The position of point P is described in 
terms of coordinate system 2 by the vector Given the vector ( X q)]^ 

from 0^ to 0^ and the angles between the axes (or lines parallel to 
them), we desire to find the vector from 0^^ to P. This vector is 
easily found by the following matrix multiplication: 



f 

4 


10 0 0 


^1 


cos (X2,Xj^) Cos (Y 2 , X^) Cos (Z 2 , X^^) 

Eq. 3 

^1 


(Yq)j^ cos (X2,Yj^) Cos (Y 2 , Cos (Z 2 , 


"l 

- ^ 


(Zq)j^ cos (X 2 , Zj^) Cos (Y 2 , Cos (Z 2 , Z^) 




The cosine terms are the cosines of the angles between intersecting lines 
parallel to the indicated axes. The sign convention is arbitrary for 
these angles since the cosine is an even function. 

We are Interested in coordinate transformations of two special types. 
One of these is the transformation due to joint angles and displacements. 
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The other transformation is due to the deflection of arm segments 
under loading. The former has been described for both rotating and 
sliding joints by J. Denavit and R. S. Hartenberg (1)* in terms 
of four independent parameters. The transformation for simple beam 
flexure, compression, and torsion will be developed in this paper. 


Transformation of Coordinates Due to Elastic Deformation 


The information we seek is the displacement and rotation of an 
arm, or more generally a jointed beam, due to the application of loads. 
The end of the beam can be described in a fixed reference coordinate 
system if one knows the transformation between the coordinate systems 
which are fixed to the individual segments. As seen in Fig. 4 the 
point p at the end of the beam can be described by two transformations, 
represented by two 4x4 matrices. The transformation relates 
system i', the end point before deflection, to system i-1. The trans- 
formation relates system i to system i'. 


Eq 4 



\ ^1 



A. E^ 


1 

0 


where: 



i-1 


the position of the origin of system i in terms 
of system i-1 


E, 


0 


transformation with no deflection 
transformation due to deflection 
a 3 X 1 vector whose elements are zero 


X 




location of point p in i coordinates = origin of i 
in this case 


* A reader consulting this paper should be aware of the fact that 

in that paper is defined with opposite sign convention of this paper 
and later papers by Denavit and Hartenberg. 

(1)J. Denavit and R. S. Hartenberg; "A Kinematic Notation for Lower- 

Pair Mechanisms Based on Matrices'* Journal of Applied Mechanics 
June 1955 pp 215-221. ^ ^ 
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Any number of these transformations may be combined by multiplying 
the transformation matrices. In terms of the reference system 0^ the 
end of a beam with N joints is located at X as is given by: 


_ 1 _ 

— A- E- An ••• A- E, A«* 

112 11 N N 

1 

- N,o_ 


0 


We would like the variation of this position vector due to applied 
forces and moments. First the elements of the E matrices must be 
found. From Eq 3 



1 


0 


0 


0 

II 

•rl 

<Vi' 

cos 

(X^.X^,) 

cos 


cos 


Eq 6 

<Vi- 

cos 

(Xj^.Y^.) 

cos 

(Yj^.Y^,) 

cos 

(Z^.Y,.) 

i 

_<Vi' 

cos 

(X^,Z^,) 

cos 

(Y^,Z^,) 

cos 

(Z^.Z^,) 


For small deflections and small angles the elements of this matrix 
simplify as follows: 



1 

0 

0 

0 


Ax 

1 

cos (90 +0 

cos (90- 0^) 

Eq 7 

AY 

cos (90- 0^) 

1 

cos (90 +0^^ 


AZ 

cos (90 + 0^) 

cos (90 ~ ® 

1 


where 0^^, 6^ and 0^ are the angles of rotation about the X, Y and Z 

axes respectively. For small angles the angles behave very nearly as 
vectors, thus the order of occurance is Irrelevant. Furthermore the 
small angle assumption allows further simplification to 
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Eq 8 


1 

Ax 

j 

AY 

AZ 


0 

1 


®z 


0 


- 8„ 


1 

e. 


0 

e 

- 0 


But these elements were expressed In terms of forces and moments 


in Eq. 1. Thus E^ may be expressed as 


Ei=| 


Eq 9 


Qt F 
"C Xii 


0 

1 




, - “bfAh- “eMl»Zlli - 


°‘xFi^Yii “xMi^Zli|“eFi^Yll'^“eMi^ii 


“Ti^Sai 


XFi^Zii ■ “xMl^iii°‘0Fi^Zii"“0MAii| 


where 

^Xii’ ^Yii’ ^zii " beam i, in terms of coordinate 

system i 

^ii’ ^ii’ ^\ii*° beam i, in terms of coordinate 

system i 


Now one must determine the forces and moments on segment i which 
result from the loads on the end of the beam. This is done in the 
following section: 

Equilibrium Forces on the Arm Segments 

k free body diagram of the beam segments between coordinate system 
i and system N is shown in Figure 5. Equilibrium requires: 
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Eq 10 

a) 2 

‘ F = 

0 

= ^Oi 

-No 

?ii 



b) 2 

M = 
, \ 

0 

= ^Oi 

”no 

"li ^Oi 

F - M 

- No -11 

where: 

lii “ 

the vector from 
terms of system 

system i 
i 

to the end 

of the arm 


= the force vector acting on the beam to the left of 
system i in Figure 5, expressed in system i 

. = the moment vector acting on the beam to the left of 

system i in Figure 5, expressed in system i 

= 3x3 rotation matrix from system 0 to system i 

Fj^q = applied force at the end of segment N, expressed in 
base coordinate frame 


^o 


applied moment expressed in the base frame 


Vectorially eq (10) may be expressed as 


Eq 11 


?ii 


^^li 


Oi 


^ii* ''oi 


Oi 


^o 


M 

No 


where r^^ x may be represented by the matrix multiplication 


-r 


Zii 


'Yii 


r xR 
ii 


Oi = 


Zll 


-r 


Yii 


-r 


Xii 


Xii 

0 


Oi 


In the above manner we can obtain the forces on the arm segments 
resulting from the loads on the end of the arm. It remains to evaluate 
the deflection of the arm by using these values in conjunction with the 
transformation matrices. 
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Arm Deflection with Load 

Having described the position of the end of the arm (after loading 
has been placed on the end of the arml^^by the coordinate transformation, 
one could subtract from this vector the vector describing the position of the arm 
before loading as in equation 13. Tlieoretlcally this would be correct. 


Eq 13 i X . [( F^) - A^ A. 




In practice the difference of these two vectors will be much smaller than 
the vectors themselves, leading to inaccuracies when the calculation is 
carried out with two few significant digits. A more practical way is to 
evaluate the partial derivative of the position of the end with respect 


to end point loads, for example and 


Eq 14a ^5)0 

5'^xno 


^2 ^2 * * * 




Eq 14b ^^0 

^^0 


i-4 


^ h '‘l 




One will now recall the assumption that the joints remain rigid. 
Because of this: 


= 3 A, 


i = 1, 2, ... , N 


^^0 


If one found that this assumption was not valid it would, be relatively 
simple to evaluate these partial derivatives and include joint flexibility. 
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By the chain rule 


~w~ 


XNO 


Ai dE^ (A2 E2 ... Ejj) 


3F 


XNO 


B-] 


^ ^1 ' 


3F. 


<A2 ^2 ... An 


XNO 


[-H 


Continuing this differentiation one eventually arrives at: (for example) 


Eq 17a ^-NO 


3F 


XNO 


S Ai • • -A-i 

i=l 


^1 ^i + 1 ^ 


3F 


XNO 


[r] 


and similarly for the other force components, as v;ell as for the 
moments: (for example) 


Eq 17b 


9^10 


N 

2 

Li =• 1 


• • *A^ 


^^i ^i+i ’Sj 


^\no 


B-] 


Then deflections are obtained as “ ^^0 ^^XNO example 


3F 


XNO 


In order to proceed we must evaluate 


3Ei 3e^ 3e^ 3e^ 3e^ and 3 e^ 

^^XNO ^^YNO ^^ZNO ^^0 ^^0 ^'^ZNO 


To do this we take the derivative of the individual elements of 
Eq . 9 as follows : 
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Eq 18 


3E. 
1 

3F 


XNO 




3F 


XNO 


-( ) 


-( ) 




“XFl ^^Yii + “xMi 


3F 


XNO 


3F 


XNO 


“XFi ^^Zii - “xMi 


3T 


XNO 


“0F1 ^^Yii + “gmI ^^Zii 


3F 


XNO 


“OFi ^^Zii 


3F 


- 0M1 


XNO 




^^XNO 1 


3F 


XNO 


3F 


XNO 




0 


-( ) 


^ ; 


I “t 


3F, 


XNO 


Similarly for F^j^q, F^jjq, ^SjNO’ ^NO’ \no* that the 

derivative of the rotation submatrix is antisymmetric 
There is but one thing left to evaluate, that is 


'!ii 


3F 


NO 




31x0 


and 


3^ii 

3^0 


Referring to Eq. 11 


it is seen that these partial derivatives are readily evaluated 
if one assumes that and r^^ x R^^ are essentially independent 
of the loading^ which they are to first order. Then 


3F 


XNO 


?ii 


\ 

Si 1 ° 

lii'^ Si 1 


Oi 


T 

0 

0 

0 

0 

. 0 . 


Eq 19 
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In general 


Eq 20 






[ 


^^0 ^^YNO ^^ZNO ^^0 ^^SrNO 


ZNo] 


L ^oi j n 

t_ -ii * ’'^Oi r ^OiJ 


These values can be substituted into Equation 18 to yield the 
derivative of the elastic deflection transformation matrix with 
respect to the end of arm loads. It has already been pointed 
out how the displacements are computed using these transformations. 
The next section will show how to arrive at the rotation of the 
end of the arm due to the loads. 


Bad Point Rotations Under Loads 

One would also like to know how the end of the arm rotates with 
applied forces and moments, i.e., determine the elements in the 6x3 
matrix C. 


' 0 * 


Eq 21 


3F 


XNO 


3F, 


YNO 


3F, 


ZNO 




®XN0 ®YN0 ®ZN0 
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Most of the work to obtain these terms has already been done. It remains 
for us to recognize the result and transform It Into the proper coordinates. 
What we have done up to now Is to take the derivative of the equation 

E, 22(a) jy - [A--] 

with respect to forces and moments applied to the arm. It was assumed In 
Eqs. 4 and 5 that the force was applied at the origin of the co- 
ordinate system, which was also the endpoint of the arm; thus ^ » 0. 

Let us use a general nonzero ^ and rewrite Eq. 22 as 

Eq 22(b) 1 0 0 0 

»0>S0 

which Is the same as 
Eq 22(c) ril Q J 

Here Is the vector from the origin of system 0 to the origin of 

system N, expressed In system 0 coordinates, while Is the vector 
from the origin of system N to the point of application of the load, 
expressed In N coordinates. 

Now we express the vector In coordinate system N' whose axes 
are parallel to the axes of system zero before loading but has the same 
origin as system N. (See Fig. 6.) The components of this new vector 3^, 
are found from the expression: 

Now 
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Since the 0 and N* axes after loading are nearly parallel, for small 
deflections cosCSO-S^^^) ° simplifications can be made as 

follows 


£q 25 


®N'0 “ 


"1 

®ZNO 


"®ZNO ®YNO 
1 “®XNO 


■®YNO ®XNO 1 


Thus to get the partial of the angles 6 above with respect to a force 
or moment , say Fjq^q one must simply multiply as follows 


Eq 26 


3F, 


XNO 


1 0 


0 

-e 


Y B 1 

NO ZNO 


ZNO ®YNO 
-0 


XNO 


^NO "^YNO ®XNO 1 


,T 1 





XNO 


1 O 


For small deflections evaluated before loading from only 

the joint angles. 

Compliance Matrix and Spring Constant Matrix 

Now we are able to piece together the above derivation to reach our 
original goal: a compliance matrix of the arm under force. Equations 

16, 17 and 26 are evaluated (as well as the similar equations for the 
other forces and moments) and one can construct the following matrix 
equation. 
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The subscripts on the matrices are understood to apply to each element. 
Due to the nature of the problem the matrix will be symmetric. The 
inverse of the matrix will be the spring constant matrix and 
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F 

c-1 

= ^NO 

AX 

o 

II 

AX 

Eq 28 

M 

NO 


NO 

A6 


CjjQ will be nonsingular for all physical cases. For some arm con- 
figurations and parameters the inverse may require excessive accuracy, 
and hence be uncalculable. In this case one must eliminate one or 
more of the directions from consideration to get an invertible matrix. 


Linear Beam Vibrations 


Up until this point we have been considering the displacements of 
and loads on a static beam. If one considers a rigid mass and inertia 
placed at the loading point of the beam, the forces and moments on that mass are 
the negative of the forces and moments on the beam. These forces and 
moments can be determined from the spring constant matrix and the 
deviation of the mass from the equilibrium position. Since structural 
damping is small, the natural frequency of the spring-mass system as 
well as the amplitude ratios of the various modes of vibration can be 
determined. Nonlinearities such as Coriolis accelerations and centripital 
accelerations can be neglected for angular velocities which are appropriately 
small. This seems to be the case in practical arm problems with small 
vibrations. The equations of. motion are then written as 


Eq 29 


where: 


M 

0 

0 

0 

0 

0 

0 

M 

0 

0 

0 

0 

0 

0 

M 

0 

0 

0 

0 

0 

0 

^XX 

'XY 

^XZ 

0 

0 

0 

^XY 


"yz 

0 

0 

0 

'xz 

^YZ 

^zz 

M = 

the 

lumped 

mass 

at the 

end of 



dt^ 


Ax 

II 

o 

Ax 

A0 


A0 

- 

NO 

- 
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^XX’ ^YY’ ^ZZ ~ mass moments of inertia of the lumped 

inertia at the end of the arm about axes parallel to the 
reference axes but through the center of mass 

Ivv ~ cross moments of inertia about axes 

AY AA Y La 

parallel to the reference axes but through the center of mass. 


For convenience Eq 29 will be rewritten as 

NO 


This can be written in state variable form as 


Eq 31 



— 





• 

d 

AX 


0 

II 

AX 


AX 

dt 









A0 

= 

K 

oj 

A0 

^ A 

A0 


AX 



AX 


AX 





% 


i 


AS 



A0 


A0 







L — J 


The dot above AX and A6 indicate a derivative with respect to time. 
The roots of the equation 


Eq 32 I si - A 1 = 0 


are the natural frequencies of the system. The amplitude ratios can 
found as for any undamped linear system. 
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Extensions - More Than One Lumped Mass 

The case of the unloaded or lightly loaded arm is one in which 
the dynamics of the arm vibration are not dominated by one lumped mass. 

The criteria for modeling with lumped masses will not be discussed 
here, but rather the use of the technique developed will be extended 
to include any number of lumped masses. Figure 7 shows schematically 
a model that one may be interested in. 

Initially one obtains spring constants between each mass point and 
its adjacent mass points. The nonequilibrium forces on each mass depend 
only on the difference in the vector positions between it and its neighbors. 
Thus for the example in Figure 7, with some change in notation: 


E, 33 J, (X^_ -X^) + ^ (X^^^ - X^) 


where 


1S.1 ■ 

L^i-lio 


= position and angular orientation for mass i, 
measured from equilibrium in base coordinates 


= spring constant matrix between mass i and mass i-1 


^ = spring constant matrix between mass i and mass i+1 


'i 

i 


= the inertia matrix for mass 1 

- 1,2,...,M 

This equation can be written for all M masses. The end masses are 
special cases 


Eq 3A 

"l5i = 

^10 -1 ^21 


Eq 35 

• • 

J X = 
M -M 

M-1 “ 

V 
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If we assemble these into one matrix expression, its form is: 


Eq 36 





J*^ K 
J 2 i<21 


I 


1- 


!-j;^(k,,+k^,)I k3^|. 

-i" 


31 2rj “2 

r 


J 1 ! [.. 

j ‘^i"\,i-i| -"i'(^i+i,i^^i.i-i^ i'^i\+i,if • 
^ 



Simplifications - Some Moments of Inertia Insignificant 


One or more of the moments of inertia of a lumped mass-inertia 
may be insignificant with respect to the mass and the other moments 
of inertia. In this case it is desirable to reduce the number of 
state variables by two by ignoring the associated angle and angular 
velocity. The moments will be continuous in the beam for the axes 
associated with the trivial moments of inertia. The other moments 
and the forces in the beam undergo a discontinuity in our lumped mass 
model due to the inertial loading. Let ^ designate the spring constant 
matrix of the entire arm, considering all points of loading . Its fotnm 
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is similar to the large matrix in Eq 36, but the J ^ terms are re- 
moved . Then 


Eq 37 


C = t 


-1 


X 

s 

F 

— 

c 

— 

£ 


M 


where X 

% 

E 

0 


unknown, possibly nonzero loading 

displacements or angles associated with the T elements 
loading terms which will be Identically zero 
angles associated with the ^ elements 


Eq 39 

X 

“ C 

F 

s 

hil 



1' 


i. 


. Si ' 

C22 


Eq 40 ^ 11 - 


Eq 41 0 = C 21 1 


'f' 

= K 

Y 

= K 

‘x 







‘ 0 ” 


0 


21 q 


\lJih2_ 


K21 I K22 



Eq 43 
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Eq 44 F = ^ \2 ^ 21 - 

Eq 45 F = (I - ^21^"^ ^1 ^ 


The above operations assume the inverse can be performed. 

The reduced equations of motion are then: 

Eq 46 M d^ Y = (I -T., C )"^ %. T 

^ — j - 12 21 11 ~ 

dt 

where M is the reduced inertia matrix obtained by eliminating the 
appropriate rows and columns from the unreduced inertia matrix. 


Example Problem 

In order to Illustrate the theory 
presented above, a computer program was developed to evaluate the 
ccmpliance matrix for an example arm. The compliance matrix was 
then input to an existing matrix manipulation program along with an 
inertia matrix to develop the equations of motion for a simple case. 

As a realistic example the arm parameters and configuration were 
taken from a proposal by the Ttartin Marrietta Company for a boom 
for the space shuttle. These are shown in Table 1. Figure 8 shows 
the arm in the configuration of the example and the distribution of 
the 65,000 lb. load. These joint angles were chosen because they 
realistically duplicate a position in a retrieve maneuver for which 
the arm might be used. It also enables a separation of modes reducing 
the number of state variables to six. This is due to the planar 
motion of the mass. Figure 9 indicates the oscillations resulting 
from an initial displacement of ten Inches in the Y direction at the 
endpoint . 
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The computer program required 0.08 hours of IBM 1130 computer time 
to evaluate the compliance matrix for six joint angle positions. This 
includes some compilation and program listing time, and the program 
could be considerably streamlined. 


Rotational Compliance at a Point 


As developed previously the transformation of coordinates due to 
deflection is given in Eq. 8. If only rotations at a point are of 
Interest the form of the transformation for small angles is : 


Eq 47 


D 


1 0 


0 1 

0 1 - 6 ^ 


Here 0^^, 0^, and 0^ are the angles of rotation due to loading about the 
X,Y, and Z axes in any coordinate system of Interest. 


These angles may be expressed in terms of the components of the 
moments acting on the point expressed in the same coordinate system as 
the angles, and a rotational joint compliance about each axis, here 
denoted and 


Eq 48 


D - 


10 0 0 

° ^ "“jz”z “jy”y 

° “jz“z ^ ”®jx“x 

° “®JY^ ^ 


This matrix can represent a joint compliance. It can then be used to 
evaluate overall arm compliance in a manner similar to the matrix E of 
Eq 8. which represents the link compliance. Note that for joint com- 
pliance due to bearing supports, etc., on the end of two adjoining links 
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which change orientation with the joint angle, two point compliance matrices 
are necessary to properly account for a change in orientation as follows 


Eq 49 


1 

= D A,D. E. 

’1 


i- i i+ i 


lr-i,i-lj 



where notation is the same as for Eq. 4 with the addition of and 

is the point compliance matrix which accounts for deflection 
of bearings, supports and drive at joint 1 which remain stationary on 
the link i-1. accounts for compliances stationary on link i. 

Flexibility and Mechanical Power Transmission 

When power is transmitted to a joint from a prime mover which is 
located away from the joint, the deflection of the link between the motor 
and the joint will depend on the manner in which the power is transmitted. 
The torque which is taken from the joint is transferred to the prime 
mover in various ways and the manner in which this is done affects the 
state of stress in the intervening segment. For example, a band drive 
with no reduction completely removes the component of moment along the 
axis of the joint and increases the compressive normal stress. A bevel 
gear and shaft drive as shown in Fig. 10 with no reduction retains the 
moment along the joint axis but shifts the moment about the link axis 
by an amount divided by the distance of transmission, A flexible 
cable drive as shown in Fig. 11 removes joint axis moments while altering 
forces mutually perpendicular to the joint and link axes and moments 
along the link axis. The effects of a particular transmission system 
must be determined by equilibrium considerations and possibly deflection 
considerations. Once determined, the effects can be represented in a 
transformation matrix which enables one to conveniently determine the 
overall compliance. For example, consider the schematic in Fig. 13. 

The segment of the arm from the fixed mounting to the motor M would be 
described by a simple beam deflection transformation matrix E of the 
form of Eq 9. The segment between the motor and gearbox would be 
described by a drive deflection transformation such as those displayed 
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in Fig. 10, 11 or 12, depending on the type of drive employed. The joint 
Itself will have a rotational compliance which accounts for the bearing 
supports, joint shafts etc. In addition servo motor compliance (although 
nonlinear) , twisting of drive shafts or stretching of drive cables, 
and deflections within the gearbox will manifest themselves in the joint 
compliance. These are rotational compliances which manifest themselves 
at a point , and should not be confused with the distributed deflection 
described by E^. For Fig. 13 the complete transformation expression for 
end point loads would be: 


"1 

Y 

^l^d°2-''^2®2+^2 




L-J 


Design Analysis and Tradeoff Studies 

Initially the simple but general case of an arm with two links and 
one joint is being studied. The criterion initially considered will be: 
maximize minimum resonant frequency and minimize the static deflection 
while penalizing the design weight. Fig. 13 shows the general case being 
studied. Even this simple case will have an unmanageable number of 
variables without certain assumptions. Among the assumptions being 
made are : 

1) Hollow circular cross sections for all arm segments 

2) Constant cross section over arm segment lengths 

3) The same homogeneous material is used in all arm segments and 
power transmission members. Nonhomogeneous materials such as filament 
reinforced composites are excluded for the time being. 

Additional assumptions will undoubtably be made as study indicates their 
reasonablemess . 

Among the questions being addressed are: 

1) What is the most desirable location of prime mover and speed 
reduction for varying sizes of prime movers? 

2) What is the most desirable allocation of structural material 
between the arm and power transmission members? 
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3) How do these decisions depend on the penalties for natural 
frequency, deflection, and weight? 

4) How do these decisions depend on the relative proportions of 
the arm? 

5) How do these decisions depend on the mode of power transmission? 

6) What Is the limiting component of the design In terms of load 
capacity? 

The analysis will be done In nondlmenslonal variables to allow the 
broadest application and the presentation of results will be graphical 
whenever possible. 

Analytical Expression for the Compliance Matrix 

An analytical expression for the compliance matrix of a two link, 
one (compliant, revolute) joint arm has been derived. Fig. 14 displays 
the case and explains the variables. Eq. 31 gives the analytical results. 

This was accomplished by using the coordinate transformation 
equation and Its derivative with respect to force. The matrix manipulations 
were carried out manually with the terms being analytic expressions 
Instead of numerical values which could be substituted In for a particular 
case. These results should avoid numerical evaluation of the compliance 
In many cases to be studied, and allow straightforward substitution of 
the arm parameters into an expression for the compliance. 

Future Work 


Preliminary work has developed the controllability matrix for the 
general case with joint angle position control. This has been used to 
show that the example problem above Is controllable using two of the 
joints. Optimal control theory can now be used to determine suitable 
feedback gains if one has access to the state variables. The state 
variables can be partially measured and partially reconstructed using 
the measured variables. Measurements might be performed via accelero- 
meters, optically, or in some other fashion. In all this future work 
the method developed here will make the determination of the equations 
of motion for arm vibration practical, even for complicated arm con- 
figurations. 
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fig 13. Scb ema'^K o-f Cose Tro(jeo^-( Stud.^ 
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